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O ' Abstract 

Motivated by recent CELO measurements and the progress of the theory of B decays, 
we investigate B PV{P = 7r,K. V = K*,p,Lo) decay modes in the framework of QCD 



improved factorization. We find that all the measured branching ratios are well accom- 
modated in the reasonable parameter space and predictions for the other decay modes 



■ are well below the experimental upper limits. We also have calculated CP asymmetries 

in these decay modes. 

PACS Numbers: 13.25Hw, 12.38Bx, 12.15Hh 
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1 Introduction 



It is of great interest and importance to investigate the decays of B mesons to charmless final 
states to study the weak interactions and CP violation. In past years, we have witnessed many 
experimental and theoretical progresses in the study of B physics with the observations of many 
B charmless decay processes and improvements of the theory of B decays. 

Most of the theoretical studies of B decays to pseudocalar and vector final states are based 
on the popular Naive Factorization approch[l|]. As it was ponited out years ago in Ref. [0], the 
dominant contribution in B decays comes from the so-called Feynman mechanism, where the 
energetic quark created in the weak decay picks up the soft spectator softly and carries nearly 
all of the final-state meson's momentum. It is also shown that Pion form factor in QCD at 
intermediate engery scale is dominated by Feynman mechanism[^, From this point, we 

can understand why the naive factorization approch have worked well for B and D decays, and 
the many existing predictions for B decays based on naive factorization and spectator ansatz 
do have taken in the dominant physics effects although there are short commings. However, 
with the many new data available from CLEO and an abundance of data to arrive within few 
years from the B factories BaBar and Belle, it is demanded highly to go beyond the naive 
factorization approach. 

Recently, Beneke et.al., have formed an interesting QCD factorization formula for B ex- 
clusive nonleptonic decays]^, |^. The factorization formula incorporates elements of the naive 
factorization approach(as leading contribution) and the hard-scattering approach(as sublead- 
ing corrections), which allows us to calculate sysmatically radiative (subleading nonfactorizable) 
corrections to naive factorization for B exclusive nonleptonic decays. An important product 
of the formula is that the strong final-state interaction phases are calculable from the first 
principle which arise from the hard-scattering kernel and hence process dependent. The strong 
phases are very important for studying CP violation in B decays. Detail proofs and arguments 
could be found in0. Here we recall briefly the essence of the QCD factorization formula as 
follows. 

The amplitude of B decays to two light mesons, say Mi and M2, is obtained through the 
hadronic matrix element {Mi{pi)M2{p2)\Oi\B[p)) , here Mi denotes the final meson that picks 
up the light spectator quark in the B meson, and M2 is the another meson which is composed of 
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the quarks produced from the weak decay point of b quark. Since the quark pair, forming M2, 
is ejected from the decay point of b quark carrying the large energy of order of nib, soft gluons 
with the momentum of order of Aqcd decouple from it at leading order of KqcD/^b in the 
heavy quark limit. As a consequence any interaction between the quarks of M2 and the quarks 
out of M2 is hard at leading power in the heavy quark expansion. On the other hand, the light 
spectator quark carries the momentum of the order of Aqcd, and is softly transfered into Mi 
unless it undergoes a hard interacton. Any soft interaction between the spectator quark and 
other constituents in B and Mi can be obsorbed into the transition formfactor of 5 ^ Mi. 
The non-factorizable contribution to -B — * M1M2 can be calculated through the diagrams in 
Fig.l. 

In this paper we study B PV decays within the framework of the QCD improved 
factorization approach |3- In Sec. II We present notations and calculations. In Sec. Ill we 
compare our numerical results with the experimental data presented by CLEO collabaration 
P]. We find that all the measured branching ratios are well accommodated in the reasonable 
parameter space and predictions for the other decay modes are well below their upper limits. 
We also give our predictions of direct CP asymmetries and time integrated CP asymmetries 
in the decay modes. Large direct CP vilation asymmetries are predicted for the decay modes 
5° TT^u, K-p+, K^uj, B^ K~K*^, K^p^, tt°K*~, K^oj. The direct CP asymmetries in 
the observed decay modes B^ 7r^p°, -n^uj and B^ n^p^ are predicted to be around few 
percentages. 



2 Calculations 



The effective Hamiltonian for B decays is given by 



/ 2 10 \ / 2 10 
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with the effective operators given by 



0]^ 


= {qaUa)v-A " {Uf3bf3)v-A, 




= (?«M^))y-A 


■ {U(3ba))v-A, 


01 






= {qoLCfi)y-A ■ 


{Cf3ba)v-A, 


Os 


= iqaba)v~A-J2q'{q'pq'i3)v-A, 


04 


= {qJ>i3))v-A 


■ Y.q'{q'pq'a)v-A, 


O5 


= {qaba)v~A-J2q'{q'i3q'/3)v+A, 




= {qabl3)v-A ■ 


T.q'{.q'l3q'a)v+A, 


O7 
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O9 
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■ Y.q' dq'iq'pq'a )v-A 
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= {gs/Sir') m, 4 a'^'' R b^ G% . 
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Here q = d,s and {q'e{u, d, s, c, b}), a and P are the SU{3) color indices and A^^, A = 1, 8 
are the Gell-Mann matrices, and denotes the gluonic field strength tensor. The Wilson 
coefficients evaluated at /i = scale are|^ 

Ci = 1.082, C2 = -0.185, 

C3 = 0.014, C4 = -0.035, 

C5 = 0.009, Ce = -0.041, 

(3) 

C7 = -0.002/137, Cs = 0.054/137, 
Cg = -1.292/137, Cio = -0.262/137, 
Cg = -0.143. 

The non-factorizable contributions to B M1M2 can be calculated through the diagrams 
in Fig.l. The results of our calculations are sumarized compactly by the following equations 



= ^ E VpqVpbl a\{BM^, M2){qu)v-A ® {ub)v-A + al{BM^, M2){uu)v-a ® {qb)v-A 

V ^ p=u,c ^ 

+ag(5Mi, M2)iq'q)v-A ® {qb)v-A + al{BMi, M2){qq)v-A ® {q'b)v-A 
+al{BMi, M2){q'q')v+A ® {qb)v-A + al{BM^, M2){-2){qq')s+P ® {q'b)s-P 

+a^r{BM^, M2)^e,,(gV)y+A ® {qb)v-A (4) 

+ (-2) (af (SMi, M2)^eg, + a^a{BMi, M2)) (gg')5+p ® {q'b)s-p 

+al{BMi, M2)^eg,{q'q')v_A ® (#)y-A 

+ {a%{BM,, M2)hq. + a?oa(5Mi, M^)) (gg')y-A ® (g'&)y-A}, 
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where the symbol ® denotes {MiM2\j2® ji\B) = (M2|j2|0)(Mi|ji|i?). Mi represents the meson 
which picks up the spectator quark through this paper. For Mi is hght vector meson and 
M2 is hght pseudoscalar meson, the effective af s which contain next-to-leading order(NLO) 
coefficients and 0{as) hard scattering corrections are found to be 
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q' = u,d,s and f = u,d, s,c,h. Cp = {N"^ ^1)/ (2A^) and = 3 is the number 
internal quark mass in the penguin diagrams enters as s j = m'j / ml . x = 1 — x 



.12 In A _ 18 + /^ + /^^, 
rrib 



(6) 



5 



fp 

Gp,g 



dxg{x)(l)p{x) 



1 - 2x, 
3 m X — 3z7r, 



1 — X 



47r^ fvfB [\ (t^Bjz) 



B 



z Jo 



X Jo 



dy- 



nl 2 

— / dx—(j)p(x), 
Jo X 

In 1-4 / dx(f>p{x) / du mm In 

3 3 mft JO JO 



rrib 
1 3 

dx-(f)p{x) 
2 



Sq — uux — %e\ 



3 

'2' 



1 i(i /"^ /"^ 

- — In h 3 / dx(j)%(x) / du nnln [s„ — nua; — ie] 

3 mfc Jo JO 



For Ml is pseudoscalar and M2 is vector, the co-efficents are 



al{BP,V) 
al{BP, V) 



0, a- (5P, y) = V),i = 3,5, 7, 8, 9, 10, 8a, 10a. 
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iV' 

iV 4^^ 
CV 
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' iV 4^^ 



Cg-FV, 
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where 



3 3 
(Cg + Cio)- 51 e/Gy(s/) + Cg-{egGvisg) + ebGv{sb)) 

f=u 



Fy = _i21n-^-18 + /f + 

nib 
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rl 1 — 2X 

fv = J dxg{x)(t)v{x), g{x) = 3 ^_ ^ Inx - 3z7r, 



N /f ^^(0)M| Jo z Jo X Jo ^ y 
Gv,g = - dx-(j)v{x), (15) 

JO X 

Gy(s„) = In 1-4 / dx(j)v{x) / mm In [s„ — tiwx — ie] . (16) 

3 3 m;, JO Jo 

Where (j)p{x) and 0p(x) are the pseudoscalar meson's twist-2 and twist-3 distribution amph- 
tudes (DA) respectively, (pvi^x) = (j)v,\\ is the leading twist DA for the longitudinally polarized 
vector meson states. We have used the fact that light vector meson is mainly longitudinally 
polarized in B decays and the tranversly polarized state is suppressed by the factor My /Mb- 
Further more m B ^ PV decays the vector meson should be completely polarized in longitu- 
dinal direction due to the angular momentum conservation requirement. In the derivation of 
the effective coefficients a^'s we have used NDR scheme and assumption of asymtotic DAs. The 
infared divergences in Fig.l.a — d are cancelled in their sum. 

With the effective coefficients in Eqs.(|]) and (|T2|) , we can write down the decay amplitudes 
of the decay modes (we only list four decay modes which have been detected in experiment 
here. The other decay modes are given in appendix A) 

MiB^ ^ ttV) = ^/pM|Fo^^-AK.fc|i?„e-'^ [ai(fi V, p-) + <(5 V, p") 



aio(5°7r+,p-) + aioa(5°7r+,p-) - a^(E%+, p") + aio(5%+, p") + aioa(5%+, p") h (17) 



MiB^' 



TT p 



V2" 



+ (<(5V, TT-) + a8(5V, vr") + a8a(5V, ^~)) + «io(5V+, ^r") + aioa(5V+, vr" 
- <(5V+, TT-) + (a^(5V, TT-) + a8(5V, TT-) + a8,(5V+, vr-)) R^^ 

+a,oiB'^p-^, 7T-) + aioa(5V, kl8) 



MiB- 



TT p 



a2{B-7i, p") - a'UB-TT, p') + - [ar{B-Tr, p'] 
+a9{B~7r, p°)) + ^aio(5~7r, p°) - a^a 
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3 1 

-al{B~7T, p°) + - (07(5-71, p°) + agiB-TT, p°)) + -aio(5-7r, p°) - a^a 

+ 9Zf^M'sA^-PXV,h[Rue-'^ [ai(5-pO,7r-) + <(fi-pO,7r-) 
+ (a6(5~P°, TT") + as{B-p\ TX-) + asa{B~p\ tt")) i?^- + aio(5-p°, vr") + aioa 
- a2(5-p°, TT-) + (a^(5-p°, TT-) + a8(fi-p°, tt-) + agal^-p", tt")) R^- 

+aio{B-p\ TT-) + aioa(5-p°, tt-)! 1;(19) 



+ (a6(5 

+03(5-0;, TT-) + a^a{B'^, TT-) j i?^- + aiQ{B~UJ, TT-) + aiOa(5~t^, TT" 

- 04(5-11;, vf-) + {al{B~uj, it') + as{B~uj, tt~) + asa{B~uj, vr-)) _R, 

+010(5-0;, Tf-) + aioa(5-o;, tt-) 



+^/^M^Fo^^^AK6^ i?ue-*^ a2(5-7r-, o.) + <(5-7r-, a;) + 2 (a3(5-7r-, o^) + a5(5-7r-, 



--[ai{B IT ,uj) + ag{B tt ,0;) — aio(5 tt ,o;) + 2aioa(5 tt ,0^^ 



04(5 TT ,0;) + 2 (^03(5 TT ,o;) + a5(5 tt ,Lu)j + \^a7{B n ,Lu) + ag{B n ,lu) 

-aio(5-7r-, o.) + 2aioaiB~ -r- , u)) /2\ 1;(20) 



Where Rh 



1-AV2 1 



A I K 



^1 and i?l 



v\- ^cb, Vud and Vus are chosen to be real and 7 



is the phase of V*^. A = |Ks| = 0.2196. = -2M^:f / {mf,{mu + m^)). 



3 Numerical calculations and discussions of results 



In the numerical calculations we use 

r(5+) = 1.65 X lO-^^s, r(5°) = 1.56 x IQ-^^s, 

Mb = 5.2792^61/, = 4.8Ge\/, = lAGeV, 

fB = 0.180GeV, U = 0.133GeV, /x = 0.158GeV, 
/x. = 0.214^6^, fp = 0.21GeV, U = 0.195GeV. 
For the chiral enhancement factors for the pseudoscalar mesons, we take 



R-n 



Rk±,o = ~l-2. 



which are consistent with the values used in [T^, |r2|. We should take care for /2^o. As pointed 
out in Ref. [^], R^^o for 7r° should be —2M^ / {mb^niu + md)) and equal to i?^± due to inclusion 
isospin breaking effects correctly. 



For the form factors, we take the results of light-cone sum rule [|1^, [14 



F^-"(0) = 0.3, F^"^^(0) = 1.13F^^"(0), 

A^^P = 0.372, A^-'^* = 0.470, 

and assume A^--(O) = 1.2*<^^(0) since we find larger Aq^'^IO) is preferred by experimental 
data. 

We take the leading-twist DA (j){x) and the twist-3 DA </)°(x) of light pseudoscalar and 
vector mesons as the asymptotic form []r3| 



'^pyi^) = 6x(l — x) 



For the B meson, the wave function is chosen as |I6| , p 

(Pb{x) = Nbx^{1 — a;)^exp 



(x) = 1. 



Ml x^ 



2ul 



(21) 



(22) 



with uob = 0.4 GeV, and Nb is the normalization constant to make Jq dx(j)B{x) = 1. (I)b{x) 
is strongly peaked around x = 0.1, which is consistent with the observation of Heavy Quark 
Effective Theory that the wave function should be peaked around Aqcd/Mb- 

We have used the unitarity of the CKM matrix V*qVub + V*qVcb + V^*Vtb = to decompose 
the amplitudes into terms containing V*gVub and V*qVcb, and 

|V;,| = A = 0.2196, \Vub/Vcb\ = 0.085 ± 0.02, 

iKfel = 0.0395 ± 0.0017, \Vud\ = 1 - AV2. 
We leave the CKM angle 7 as a free parameter. 

The branching ratios of two body B decays is given by 

tb 



(23) 



Br{B M1M2] 



■\M{B M1M2] 



(24) 



IG-kMb 

For the case that the final state / is non-CP-eigenstate, the direct CP asymmetry parameter 
is defined as 



A dir 
^CP 



\M{B+^f)\^ + \MiB-^f)\^' 
9 



(25) 



and 

.... ^ \M{B' ^ fr - IMjB^ ^ f)\' 

CP |_M(50^/)|2 + |_^(^0^/)|2- ^ 

For the neutral B decaying into CP eigenstate /, i.e., / = /, the effects of — B'^ mixing 
should be taken into account. Thus the CP asymmetry is time dependent, which is given by pi 

1 + 



Acpit) = A'^i'pCosiAmt) - ^sm(Amt), (27) 



and the time- integrated CP asymmetry is obtained through 

1 - - 2/me(Am/r) 

(l + |e|2)[l + (Am/r)2]' ^^""^ 
where Am is the mass difference of the two mass eigenstates of neutral B mesons, and is 
the direct CP asymmetry defined in eg. (p6D . The parameter ^ is given by 

^ v;Mf\Heff\BO) 

^ Vt,V:Af\Heff\BO)- ^ ' 

The numerical results of the branching ratios B PV are shown in Fig. 2 as the function 
of CKM angle 7. We can see from Fig.2.1, 2.2, and 2.3 that for the three detected channels the 
predicted branching ratios agree well with the CLEO experiment data P]. Our predictions for 
other decay modes are well below their 90%C.L upper limits. 

There are several works available with deatil analyses of the CLEO new data of the decays 



of B to charmless PV states[|ri|, |T2|, |T9[. Cheng and Yang have renewed their predictons with 



the "generalized factorization" framework||ll||. It is worth to note that the shortcommings in 
the "generalized factorization" are resolved in the framework employed in this paper. Non- 
factorizable effects are calculated in a rigorous way here instead of being parameterized by 
effective color number. Since the hard scattering kernals are convoluted with the light cone 
DAs of the mesons, gluon virtuality k"^ = xml in the penguin diagram Fig.l.e has well defined 
meaning and leaves no ambiguity as to the value of k"^, which has usually been treated as a 
free phenomenological parameter in the estmations of the strong phase generated though the 
BSS mechanism I^Oj. So that CP asymmetries are predicted soundly in this paper. We present 
the numerical result of the branching ratios of 5 — PV decays in Table. ^ with the relevant 
strong phases shown explicitly. It shows that the strong phases are generally mode dependent. 

Hou, Smith and Wiirthwein have performed a model dependent fit using the recent CLEO 
data and found 7 = 114j;2i degree. Using SU(3) flavor symmetry, Gronau and Rosner have 
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Table 1: Strong phases in the branching ratios (in units of 10 ^) for the charmless decays 
modes studied by CLEO. 7 = ArgV*^,. 



B{B- 7r-p°) = 6.65|0.11e-^s6-5° + e"^^ 


2 B{B° TT+p-) = 19.79|0.11e*9-°2° + e-^T| 


2 


B{B^ n~p+) = 13.43|0.03e^i^2° ^ 


2 ^ TT-cj) = 10.59|0.065e*26-oi° + e~'^ 


2 


B{B^ _ vrV°) = 0.11|0.21e2-90° + g-^^p 


B{B- Tx^p-) = 10.81|0.176e^^-20° + e-'^' 


2 


B{B^ n°Lu) = 1.49 x 10-'^\l.QAe'^'^^° + e" 


B{B- K~p°) = 0.55|0.24e-^i62° ^ g-n 


2 


B{B- Ti-R*^) = 0.0012|56.4e-*i5-^° + e" 


5(5- K-K*^) = 0.030|2.86e*i^^° + e"^^ 


2 


B{B~ Tr°K*~) = 0.59|2.80e-^^''9° + e''"' 


|2 ^ if-cu) = 0.80|0.48e-*''-23° + e"*^] 


2 


B{B^ K^uj) = 0.72|0.81e-*"-«° + e-'^\ 


2 ^ i^"p+) = 0.96|0.63e-^^-20° + 6"^^ 


2 


B{B^ 7rOJ?*0) = .004|12.89e*67-6i° + g-n 


2 



analysed the decays of B to charmless PV final states extensivly and found several processes 
are consistent with cos 7 < 0. In this paper we find cos 7 < is favored by the B^ 7r"p° and 
B^ 7r^p+ + 7r+p~ if their experimental center values are taken seriously. To meet its center 
value with C0S7 < 0, B^ tt^uj would indicate larger form factor i.e. Aq^'^{0) > Aq^^{0). 
In our numerical calculation, we have taken A^^'^{0) = 0.446 which is still consistent with the 
LCSR results 0.372 ± 0.074|jl3l. It is also interesting to note that B^ tx'^ p~ is suppressed 
by cos 7 < while B^ 7r~p"'" is enchanced. The defference between Br{B^ ^^P~) and 
Br{B^ — i> TT^p'^) is much more sensitive to 7 than their sum. 

For comparison with the results in the literature, we table our predictions made for 7 = 100° 
in Table.|^. We find that most of our predictions agree with Ref. [Tl[]. For Br{B^ tt^K*^), 



our prediction is much smaller than the prediction of Ref. which exceeds the upper limit 
slightly. 

The direct and time-integrated CP asymmetries are shown in Fig3. and Fig. 4 respectively. 
For 7 around 100°, the direct CP asymmetries in the decay modes B~ K~K*^, K~p^, tt^K*~ , 
K~uj and B^ K"p~^, K^uo are as large as (±)10% ~ (±)20%. Unfortunatly the decay rates 
of these decay modes are quite small. From Fig. 4, we can see that the time integrated CP 
asymmetries in B^ P^tt", p°7r° , cutt" are about ±30% ~ ±60% for 7 around 100°. 
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Table 2: Branching ratios (in units of 10~^) for the charmless decays modes studied by CLEO. 
Experimential results are taken from Our results are made for 7 = 100°. Cheng and Yang's 
preferred predict ions |]11| (the case Nlf^{LL) = 2 and N^^^{LR) = 6) are taken for comparison. 
The for m factors used in Ref . [p!I|] are very similar to ours. 



Decay modes Our results Ref|[TT|| CLEO B 



Hit 



B or 90% 5 ULl 



£J 


TX P 


o.yu 
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1 n /L+3-3 _(_ 9 1 


B" - 
B^- 


4- — 


18.6 
13.5 


18.2 
14.2 


bum up 
27.61?;^ ±4.2 


bum up 
27.Qtji ± 4.2 


B- 


— > TT^Ll! 


9.82 


10.7 


11.3t|^ ± 1.4 


11.3l|^ ± 1.4 






0.11 


0.75 


1.6l?:^±0.8 


<5.5 


B- - 




10.0 


13.1 




<43 


B^ 




0.004 


0.02 


n 0+1.9+1.0 

•"••O-0. 8-0.8 


<5.5 


B~ - 


K'p' 


0.71 


1.10 


8.4l^:°±1.8 


<17 


B- - 


^ 7r-i?*o 


3.85 


3.64 


T.Qt^l ± 1.6 


<16 


B- - 


> K-K*^ 


0.26 


0.39 


n r,+i. 3+0.6 

U-U_o.O-0.0 


<5.3 


B~ - 


^ TT^K*- 


6.42 


4.34 


7.11^^1^ ±1.0 


<31 


B- - 


K~iu 


0.97 


2.24 


3.2l?:^±0.8 


<7.9 


B'- 


K^uj 


0.40 


1.89 


lo.oil:^ ± 1-4 


<21 




^ K^p+ 


1.29 


3.49 


16.0i^i±2.8 


<32 


B'- 




0.58 


3.92 


n n+i-3+0.5 
u-u_o. 0-0.6 


<3.6 
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4 Summary 



In this paper we have calculated the branching ratios and CP asymmetries of the charmless 
decays B PV{P = {7r,K),V = {p,uj,K*)) in the QCD improved factorization approach 
which have been formed recently by Beneke et. al. P, ^. 

We have used LCSR form factors F^^'^'^(O) and Aq^ (0) as inputs. The results of 
Br{B~ — s> TT^p^) and Br{B^ n^p^) agree with CLEO data]^ very well and favor cos 7 < 
if their experimental center values are taken seriously. To meet its experimental center value 
and C0S7 < 0, the decay B'^ ti~uj will prefer larger form factor ^^^'^(0). For the other 
decay modes, the branching ratios are predicted well below their 90%C.L upper limits given 
in Ref^]. We have also compared our results with Cheng and Yang's renewed results of the 
branchig ratios. For many decay modes, our results agree with theirs. 

Working in the QCD improved factorization approach, we are allowed to calculate the strong 
phases to make predictions of CP asymetries for the decay modes more sound than before. We 
find direct CP asymmetries in the observed decay modes are around few percentages level. 
Direct and time- integrated CP asymmetries in those decay modes have been tabled in Fig. 3 
and Fig. 4 respectively. 
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Appendix A 



The decay amplitudes of some of tlie B PV decays in terms of the effective coefficients 



tti s: 



3 



/.M^Fo^^-AK.^ Rue-'' a2{B\\ p°) - <(5°7rO, 



2 ^a,{B\', pO) + a,{B\\ pO)) + -a,,{BW p") - a^oa{B\\ p') 



-aliB'n', pO) + ^ (a,(5°7r°, p°) + a,{B%\ p°)) + ^aio(5°vr°, p°) - aioa(5°vr°, p° 



2\/2 



+-aio(fiV,vr°)-aioa(fiV,^' 
-| (a,(5V, tt") - a9(5V, vr°)) + ^a^o (^V, 7r°) - a^oal^V, tt")] |; (30) 



M{B- ^ ttV) = ^fpMlF,^--\V,k[Rue-'' [ai(5-7r°, p") + <(5-7r°, p" 

+aio(fi~7r°, p^) + aioa(5"vr°, p~ 
a2(5-7r°, p-) + aio(5-7r°, p") + aioa(5-7r°, p") 



(ay(5V-, 7r°) - a,{B'p~, 7r°)) + ^aio(5V-, vr°) - aioa(5V, vr°)] };(31) 
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-2 (a3(5°7r°, 00) + a5(5°7r°, 00)) + ^ (a7(5°7r°, cu) + ag{B%\ u 

-aio(5V,Lj) + 2aioa(5V,o^ 



+ ^ (a7(5°7r°, + a^{B\^, u) - ai^{B\\ u) + 2aioa) } 



2 

Gf 



(32) 



G, 



fKMlA^^PV,,{l - O.SA^)^ i?,e-*^ a^{B- p\ K') + al{B- p\ K 



+rki {al{B-p\ K-) + as{B- p\ !<-) + a8,(5-p°, A")) + aio(5-p°, iT") + ai^a{B- p\ K' 



+ 



al{B-p\ K-) + rki (al{B- p\ + asiB^p'', K^) + a8a(5"p°, /<- 

+a^,{B-p\K-) + a^oa{B-p\K-\ 



+-a,{B-K-,p') 



3 

+ 2 



a2(i?-A^-,p°) + -a7(i?-A-,pO) 
a7(5-i^-,pO) +a9(5-ir-,p°)l j;(33) 



<(S-7r-,i^*)--aio(S-7r-,i^* 



-aiOa(5 IT ,K* 



+ 



al{B Tx ,K*)--aio{B vr , K*) + aioa{B tx ,K 



;(34) 



G 



-47 



al{B-K-,K*) - -aio{B-K-,K*) 



al{B-K-, K*) - -aw{B-K-, K*) + a^^a^B' R- , K* 



(35) 
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TT 



+ 



+aio(5-7r°, K*-) + aioa(5-7r°, i^*" 
a^(5-7rO, Ji *-) + aio(5-7rO, i^*") + aioa(5-7rO, K*~) 



2 



-ag{B-K*- 



■17 



TT 



a7(5-/r-,7r°) 



TT" 



3 r 
2 



a7(5-A'*-,7r°) -a9(5-K*-,7r°)l l;(36) 



+ 



a2{B-K-, uj) + 2 (a3{B-K-,uj 



+a^{B-K-,uo)) + - (07(5-/^-, c^) + ag{B-K-,uj] 



2 (^a3(5-A'-,cj) + a5(5-i^-,cu)j + - i^ajiB" K~ ,uj) + ag{B~ ,uj] 
+Rk- {aeiB-uj, + a^{B~uo, K') + a8a(5~cj, AT^)) + aio(B"cj, iT^) + aioa(5^cu, K~ 



+ 



i^-) + Rk- {al{B-u, R-) + 03(5-0;, K-) + a8a(5-cu, i^" 

+aio(5-w, K^) + aioa(5-u;, A^" [>;(37) 



+ar>{Bm\ uo)) + - {a,{B^K\ 00) + a9(S°i?°, 



2 (03(5° J?o, cu) + a^iBm"^, iu)) + - (07(5° A^°, cu) + 09(5°/?°, c^; 



0, , i^o^ 



--a^o{B'u;,K") + a,oaiB''uj,K") + 



;(3^ 
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+Rk- K-) + a8(5V+, K 



+Rk- {a'',{B'p-^,K^ 



+aio(5V,A-) +aio„(5V,A-)J + 
a8(5V+, K~) + a8a(5V, K-)) + aio(5V+, A'~) + aioa(^V, A'^) h(39) 



-17 



a2(SOi?*o,7rO)-^(a7(50i?*°,7rO) 
3 r 



a7{B^K*'',Tr'')-ag{B^K 



Or>*0 ^On 



+awa{B'7r',K*^' 



al{B 



1 



a^[n IT ,K ) 



aio(i?V,/?*0) + aioa(i?"vr",A' 



0^0 r:'*0\ 
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Figure captions 



Fig.l Order as non-factorizable contributions in S — > M1M2 decays 

Fig. 2 BR{B — > PV) as a function of 7 are shown as curves in units of 10~^. The BR measured by 
CLEO Collaboration are shown by horizontal solid lines. The thicker solid lines are its center values, 
thin lines are its error bars or the upper limit. 

Fig. 3 Direct CP asymmetry of i? — > PV as a function of 7 

Fig. 4 Time-integrated CP asymmetry as a function of 7 with the mixing parameter Am/F = 0.723 
in the SM 
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.2.1, Br{B° K-p+ + K+p-) vs 7 Fig. 2. 2, Br{B- k' vs 7 Fig. 2. 3, Br{B- k-lo) vs 7 
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Fig. 2.4, BriB"^ ttV) vs 7 Fig. 2. 5, Br{B- n'^p-) vs 7 Fig. 2. 6, 5r(5'^ tt'^cj) vs 7 




50 100 150 200 250 300 350 




0.8 
0.7 
. 6 
0.5 
0.4 



50 100 150 200 250 300 350 




50 100 150 200 250 300 350 



Fig. 2. 7, Br{B- K*^) vs 7 Fig. 2. 8, Br{B- K' K*^) vs 7 Fig. 2. 9, Br{B- K~p°) vs 7 
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Fig. 2. 10, Br{B° K-p+) vs 7 Fig. 2. 11, Br{B° 7r°/i'*°) vs 7 Fig. 2. 12, Br{B- tt^K*-) vs 7 
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Fig. 2. 13, Br{B- R-lo) vs 7 Fig. 2. 14, Br{B^ K^lo) vs 7 
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Fig.3.1a, Ai!;{B° p+) vs 7 Fig.3.1b, Ai!;{B° n+p-) vs 7 Fig.3.2, Aj!; 7rV°) vs 7 
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Fig.3.3, Afj{B- ttV) vs 7 Fig.3.4, Afj{B- vs 7 Fig.3.5, Afj{B- 71'^) vs 7 
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Fig.3.6, A'll'iB^ 7r°cj) vs 7 Fig.3.7, A'^l' {B' tt' K*^) vs 7 Fig.3.8, A'^l' {B' K' K 
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Fig.3.9, A'll'iB- K-p^) vs 7 Fig.3.10, A'^l' {B^ K' p+) vs 7 Fig.3.11, A'^l' {B^ tt^K 
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Fig. 3. 12, Ajj''{B- TT^K*-) vs 7 Fig. 3. 13, Ajj'' {B' R-cj) vs 7 Fig. 3. 14, A^l'' {B° K^cj vs 7) 
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Figure 4: 
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